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^ ■ Abstract 

The decoherent histories approach to quantum theory is applied to a class of reparametrization 

(N . 

invariant models, which includes systems described by the Klein-Gordon equation, and by a min- 

(N 

isuperspace Wheeler-DeWitt equation. A key step in this approach is the construction of class 

m ; 

operators characterizing the questions of physical interest, such as the probability of the system 
entering a given region of configuration space without regard to time. In non-relativistic quantum 
mechanics these class operators are given by time-ordered products of projection operators. But 



o 

ON 

o 
in 
o 



in reparametrization invariant models, where there is no time, the construction of the class op- 
erators is more complicated, the main difficulty being to find operators which commute with the 
Hamiltonian constraint (and so respect the invariance of the theory). Here, inspired by classical 
considerations, we put forward a proposal for the construction of such class operators for a class 
of reparametrization-invariant systems. They consist of continuous infinite temporal products of 
Heisenberg picture projection operators. We investigate the consequences of this proposal in a 
number of simple models and also compare with the evolving constants method. 

PACS numbers: 04.60.-m, 04.60. Gw, 04.60.Kz, 03.65.Yz 
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I. INTRODUCTION 



A. Opening Remarks 



A problem attracting some interest in recent years concerns the quantization of simple 
cosmological models which possess no intrinsic time parameter, and which are described by 
an equation of the Wheeler-DeWitt type 



HV = 



The absence of a time parameter together with the associated reparametrization invariance 
represent a particular challenge to conventional methods of quantization and interpretation 
and it has proved surprisingly difficult to extract probabilities from the wave function. Two 
particular approaches have made interesting progress in this area: the evolving constants 
method , and the decoherent histories approach 

0, S 0, Q Q, Q- 

The aim of this paper is to develop further the decoherent histories quantization of these 
"timeless" theories described by an equation of the form Eq. (|l.ljl . 



B. The Decoherent Histories Approach 



We first briefly review the decoherent histories approach in non-relativistic quantum 
theory described by a Schrodinger equation jl^, (ill 12, 13 Q, Q In the decoherent 
histories approach probabilities are assigned to histories via the formula, 



p(ai, a 2 , ■ ■ ■ ) = Tr (C„pC^) 



where Ca denotes a time-ordered string of projectors at times ti ■ ■ ■ t r 



Ca — Pa n (tn) ' ' ' Pa 2 (^l) 



;i.2) 



;i.3) 



and a denotes the string a±, ■ ■ • a n . The projections operators are in the Heisenberg 
picture, 

P ak (t k ) = e lHt "P ak e- tHtk (1.4) 

where the projectors satisfy 



and 

P a Pf3 = SapP a (1.6) 

We are interested in sets of histories which satisfy the condition of decoherence, which is 
that decoherence functional 

D(a,g!)=Tr(CnpCl?) (1.7) 

is zero when a 7^ a'. Decoherence implies the weaker condition that R,eD(a, a') = for 
a 7^ a', and this is equivalent to the requirement that the above probabilities satisfy the 
probability sum rules. We normally work with the stronger condition of decoherence, which 
is related to the existence of records, corresponding to generalized measurements jl5l l21|. 

Now some simple observations relevant to what follows. The class operators Eq. (|l.HJ) 
defined above satisfy 

( L8 ) 

a 

However, in non-relativistic quantum mechanics, one can equally well define the class oper- 
ators to be 

n - P p -iH(t n -t n -i) p >t -iH(t2-ti)p (1 q\ 

as long as the initial density matrix is redefined to absorb a unitary evolution factor (the 
unitary factors at the final time cancel out in the decoherence functional). This alternative 
class operator satisfies 

^C7„ = e- iH(tn - tl) (1.10) 

a 

The distinction between these two class operators is trivial in non-relativistic quantum me- 
chanics but not so in reparametrization-invariant theories where one has to ask afresh what 
a class operator actually is. The definition Eq. (jl.9|) with the property Eq. (|1.10|) views the 
class operator as the decomposition of a propagator, and is best thought of in terms of 
a restricted sum over paths in a path integral. The definition Eq. (|1.3|) with the property 
Eq. (jl.8|) . on the other hand, views a class operator as the generalization of a projection 
operator since clearly it would be a projection operator if all the projections at different 
times commute. The difference between these two views is irrelevant in non-relativistic 
quantum mechanics but can have a significant influence when it comes to generalizations to 
reparametrization-invariant theories, as we shall see shortly. 
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C. Decoherent Histories for Systems without Time 



The structure of the decoherent histories approach is very general and readily applies to 
a wide variety of situations, provided one specifies a number of things in the construction of 
the decoherence functional Eq. (jl.7j) . such as the inner product structure and the form of the 
class operators. Here we are concerned with reparametrization-invarariant theories which 
are characterized by a constraint equation of the form 

H\V)=0 (1.11) 

where H is usually quadratic in all the momenta. Important examples are the Klein- 
Gordon equation of relativistic quantum mechanics and the Wheeler-DeWitt equation of 
quantum cosmology. The solutions to this equation are usually not normalizable in the 
usual Schrodinger inner product, so we use instead the so-called induced (or Rieffel) inner 
product [22|. This involves first considering eigenstates of H, 

H\^ Ek ) = E\^ Ek ) (1.12) 

where A; is a degeneracy label. The spectrum of H is typically continuous and in the usual 
inner product we have 

(^Ekl^E'k') = S(E - E')5(k - k') (1.13) 

The induced inner product between the eigenstates with the same E (including E — 0) is 
then defined, loosely speaking, by dropping the 5 function in E, that is 

(*Ek\*Ek>)l = S(k-k') (1.14) 

In practical terms, this means working with the usual inner product, regularizing all expres- 
sions by working with eigenstates of H and then dropping 5{E — E') at the end. 
It is also useful to note that eigenstates of H may be written 

\* E ) = 5(H-E)\<f>) (1.15) 

for some fiducial state \<p). Then 

(VeI^E') = 5(E - E')((f)\5(H - E)\(f>) (1.16) 
so the induced inner product between two eigenstates with the same eigenvalue is 

(*e\*e)i = W(H ~ E)\cj>) = (y E \<j>) (1.17) 
4 



See Refs. |9|, |10|, |ll( for applications similar to those considered here. 

We are interested in the construction of the class operators for systems of this type. The 
key property of reparametrization-invariant theories is that they generally do not possess a 
variable to play the role of time, hence all questions that one asks about the system must not 
refer to time in any way. We will concentrate on the following useful question: given that 
the system's state satisfies the constraint equation, what is the probability of finding the 
system in a region A of configuration space, without regard to time? The question is clearly 
a sensible one classically, since the system has a number of classical trajectories and one 
can ask what proportion of them pass through the region in question. Moreover, classically, 
it is also a reparametrization-invariant question, since an entire classical trajectory is a 
reparametrization-invariant object [llj|. To answer this question in the quantum case we 
need to find a suitable class operator. 

It is generally held that for reparametrization-invariant theories, the most significant class 
of physical questions involve operators which commute with the Hamiltonian Q, 0, H 2^, 



2J, |2i|. These are referred to as "observables" and are the analogues of gauge-invariant 



quantities in gauge theories. This issue is not without debate and subtlety in the case of 



theories invariant under reparametrizations 



but in this paper we will go along with this 



general idea. We therefore seek a class operator Ca satisfying 

[H,C A ] = (1.18) 

and which corresponds to the statement that the system passes through the region A without 
reference to time. With reference to the discussion of Section IB, we adopt the projection 
operator view of the class operator, so the Ca becomes the identity when A becomes the 
entire configuration space and is also a projector when everything commutes. 

On the other hand, if we looked instead for a class operator which is the analogue of the 
propagator form Eq. (jl.9j) . then Ca become 5(H) when A becomes the entire configuration 
space (although this is essentially the identity when operating on solutions to the constraint 
equation). We then expect that class operator to satisfy the constraint equation 

HC A = (1.19) 

However, this propagator viewpoint naturally leads to a path integral construction which, 
in earlier works, was found difficult to reconcile with the constraint equation, Eq. (jl.l9j) 
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lOL lllj . In this paper we will therefore concentrate on the projection operator form, the 
generalization of Eq. (jl.3|) . and reparametrization invariance is easily maintained. 

On a more practical note, in order to evaluate the decoherence functional and probabil- 
ities, we will need to evaluate expressions of the form where A commutes with H 
and \ty) is an eigenstate of H. Using Eq. (jl.l5| ). we have 

(*e\A\* e >) = W( H ~ E)A5(H - E')\cf>) 

= 5(E - E')((j)\8(H - E)A\4>) (1.20) 

so in the induced inner product we have 

(*e\A\V e )i = (<f>\S(H - E)A\(j)) 

= {^E\A\(j)) (1.21) 

For example, using this formula, the decoherence functional (with class operators commuting 
with H) is conveniently written 

D(a,a') = (*E\Ct,C a \*E)i 

= (4>\Cl5{H - E)C a \(j>) (1.22) 



D. This Paper 

In Section 2 we describe the earlier attempts to construct class operators for timeless 
models and the associated mathematical machinery that we will need here. In Section 3 
we describe our proposal for new class operators which are compatible with the constraint. 
They consist of infinite products in time of projection operators in the Heisenberg picture. 
We compare the decoherent histories approach with the "evolving constants" method in 
Section 4. In Section 5, we show that the decoherent histories approach together with the 
new class operators gives sensible and expected results for the non-relativistic particle in 
parametrized form. In Section 6 we compute the class operators and decoherence functional 
for a simple one-dimensional example and we apply this understanding to the relativistic 
particle in Section 7. We look at some simple examples in two dimensions in Section 8. The 
very different and simpler case of systems of harmonic oscillators is covered in Section 9. 
We summarize and conclude in Section 10. 
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II. BACKGROUND 



We begin by describing the propagator viewpoint for the construction of class operators 
for timeless systems. As stated, this has difficulties in relation to the constraint, but the 
details of the construction are important. Recall that we are interested in the question, given 
that the system is in an energy eigenstate, what is the probability of finding the particle in 
a region A of configuration space, without regard to time? 

We will consider a system whose d- dimensional configuration space is M d and it will 
generally be useful to denote their coordinates by a vector x, although when talking about 
the relativistic particle, we will use the usual notations x or x M . The propagator approach to 

is in the configuration 



In this 



defining the class operators is to define them by summing over all pat 
space between given end points which pass through the region A 
approach, the class operator of interest is therefore given by 

/oo 
dT A (x",T|x',O) (2.1) 
-oo 

The integrand is given by a standard path integral (of non-relativistic type) 

2 A (x",T|x',0) = j Thcexp(iS[x.(t)]) (2.2) 

where the sum is over all paths from x' to x" in time T which pass through A and ,f>[x(t)] 
is an action of the usual form 

Slx(t)]= / dtifijXiXj-Vix)) (2.3) 
Jo 

for some metric on the configuration space fij (whose explicit form will be unimportant 
in this section). This definition seems reasonable since it is an obvious generalization of 
Eqs. fjl.9j) . (jl.lOjl . Also, if we let the region A become the whole configuration space, then 
g is a solution to the Schrodinger equation, and, since the integration range of T is infinite, 
C(x",x') is a solution to the constraint equation, 

EC = (2.4) 

(If we are interested in an eigenstate of energy E, then we may assume that E has been 
absorbed into the potential V). However, there is a fundamental problem with this con- 
struction, which is that the class operator does not appear to satisfy the constraint equation 
everywhere, except for the case when A is the whole configuration space. 



To see this, it is necessary to go into more detail about the construction of the above 
class operators. These details will also be important for the construction of projector- type 
class operators which commute with the constraint. We first introduce the (provisional) 
class operator for not entering the region A, which is given by a restricted sum over paths 
that do not enter A, 

/oo 
dT # r (x", T|x', 0) (2.5) 
-oo 

Here g r (x", T|x', 0) is the non-relativistic restricted propagator, defined by a sum over paths 
in fixed time T that do not enter A. It vanishes when either end-point is in A or on its 
boundary. Since the set of all paths between the fixed end-points either pass through A or 
not, we have 

«?(x", T|x', 0) = <? r (x", T|x', 0) + <? A (x", T|x', 0) (2.6) 

and correspondingly 

C(x", x') = CV(x", x') + C A (x", x') (2.7) 
Here, C(x",x') denotes the sum over all paths, and in fact 

C(x",x') = (x'>(#)|x'> (2.8) 

There is a way of writing the restricted propagator which will be useful for later sections. 
We introduce the projection operator P onto the region A, 

P = J d d x\x.)(x\ (2.9) 

together with the complementary projector P = 1 — P onto the region A outside A. Suppose 
we divide the time interval [t',t"] into discrete points, t' = t < t± < t 2 < ■ ■ -£ n -i < t n = t", 
where tk+i—tk = St. We introduce the operator version g r (t", t') of the restricted propagator, 

r (x",t"|x , ,f), so 

g r (x",t'y,t') = (x"\g r (t",t')\x') (2.10) 
The operator version is then given by 

9r (t",t') = lim p e -imn-t n -l) . ..p e -iH(ti-to)p (2.11) 

where the limit is 5t — > 0, n — > oo with n8t = (t" — t') held constant. From this one can 
clearly see that g r (x." ,t"\x.' ,t') vanishes if either end-point is in A. One may also see that it 
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does not quite satisfy the Schrodinger equation, but satisfies instead, 

(i-^-H^j g r (t",t') = [P,H)g r (t",t') (2.12) 

Because P is a projection operator, in the x representation the right-hand side consists 
only of 5-functions on the boundary of A. So the restricted propagator almost satisfies 
the Schrodinger equation, but just fails at the boundary. Correspondingly, when used to 
construct the class operator C r (x", x') in Eq. (|2.5jl . it fails to satisfy the constraint because of 
5-functions on the boundary. Consequently, Ca also fails to satisfy the constraint, because 
of Ea. (l2~Tj) . 

It should be noted that constructions such as Eq. (j2.5J) can, in fact, be argued to be 
reparametrization invariant and one might therefore expect that it satisfies the constraint 
equation. The fact that is does not quite satisfy the constraint is related to subtle differences 
between the way reparametrizations act in configuration space versus phase space 2f| . 

Another useful formula for the construction of these class operators is the so-called path 



decomposition expansion 



27, 



28 



29j. The propagator (7a( x " ; £"I x 'j t') is given by a sum over 
paths which enter the region A. These paths may be partitioned according to the time 
t c and place x c at which they cross the boundary £ of A for the first time. The crossing 
propagator may then be written, 

<?a(x", tV, = I dt c I d d ~ l x c g{x", t"|x c , Q ■ V^(x c , £ c |x', t') (2.13) 

J t J s 

where the normal n points towards the restricted propagation region. Although the restricted 
propagator g r vanishes on S, its normal derivative does not (if defined by first taking the 
derivative and then letting x c approach the surface from within the restricted propagation 
region). In fact the combination (i/2m)n- Vg r (x c , t c |x', t') represents a sum over paths which 
do not cross £ but end on it. 

The path decomposition expansion was used in Ref.fiol] to compute class operators cor- 
responding to crossings of a spacelike surface in relativistic quantum mechanics. As stated 
above, the class operators constructed using the above methods failed to satisfy the con- 
straint. However, following a suggestion in Ref. jjj, it was shown that operators satisfying 
the constraint and yielding sensible results could be obtained by some simple and physically 
reasonable modifications. But this procedure is rather ad hoc and it was not clear how to 
turn it into a general definition of the class operator. 



Another clue as to how class operators should be constructed was found in Ref.jll|, which 
considered general minisuperspace models and attempted to construct class operators for 
them. The starting point was the construction of probabilities for timeless coarse grainings 
in the classical theory. Suppose we have a classical theory described by a phase space 
probability distribution function u>(p, x) satisfying 

{H,w} = (2.14) 

(the classical analogue of the constraint equation). Let /a(x) denote the characteristic 
function of the region A, so is 1 or 0, depending on whether x is inside or outside A. Now 
introduce the classical solution x c/ (t) passing through the phase space point p, x. Then the 
quantity 

/>oo 

ta= dt f A (x d (t)) (2.15) 



is the amount of parameter time spent by the trajectory in the region A. This quantity has 
the important property that it has vanishing Poisson bracket with the Hamiltonian, 

{H,r A } = (2.16) 

so is a classical observable. To determine whether or not the trajectory passes through A, 
we only need to know if ta is positive or zero. It follows that the probability of entering A 
is given by 

Pa = j d d pd d x w(p,x) 6 (J dt / A (x d (t)) - e^j (2.17) 

Here, e is a small parameter which goes to zero through positive values, and is included 
to avoid ambiguities in the ^-function at zero argument. The whole expression is invariant 
under reparametrizations, since each part of it is. Similarly, the probability for not entering 
the region is obtained by flipping the sign of the argument in the ^-function. (Note that 
there is an issue of normalization of w(p, x) in Eq. (j2.17|) . since w(p, x) is constant along the 
classical trajectories. This issue is in fact resolved by the normalization in the analogous 



quantum discussed in Ref.jHj). 



Inspired by the classical case, it was suggested in Ref. [ll] that in the quantum case, the 
class operator in the semiclassical approximation is given by 

Ca^Xq) = 9 (y°° dt f A (4(t)) ~ 5(x/,x ) e^ x " x °) (2.18) 
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where Be tA is the usual unrestricted semiclassical propagator, and Xg(i) denotes the clas- 
sical path connecting xo to x/. This object satisfies the constraint in the semiclassical 
approximation and gave sensible results, but no fully quantum version was given. 

Note that the classical and semiclassical results Eqs. (j2.17j) . (J2.18|) involve entire classical 
trajectories, not trajectories of finite length between fixed end-points as indicated by con- 
structions such as Eq. ()2.1|) . This is significant since, as argued previously, a whole classical 
trajectory is reparametrization invariant, whereas a section of classical trajectory is not jll|. 
Hence one of the key ideas in the quantum theory is to get away from propagation between 
fixed end-points and towards objects which capture the idea of an entire trajectory, as we 
will see shortly. 

III. NEW CLASS OPERATORS 

Given the above background and difficulties, the question now is how to define class 
operators that commute with the constraint and that give sensible semiclassical results. In 
this section we will focus on the case in which the unphysical parameter time t takes an 
infinite range. The special case in which the parameter time is periodic (bound systems) 
will be treated in Section 9. 

A useful hint towards constructing class operators comes from the ^-function used in the 
expressions (j2.17j) and (j2.18j) . Suppose we are interested in the probability of not entering 
A. Let is write the appropriate ^-function in terms of its Fourier transform, 



where ta is given by Eq. (j2.15)l . As stated, this object is reparametrization- invariant in that 
it has vanishing Poisson bracket with H . Now consider a discretized version of the time 
integral for ta, so we split the time into small intervals of size 5t, and we have 



with exact agreement with the original expression in the limit 5t — > 0. This then means that 
the ^-function is given by the continuum limit of the expression 




(3.1) 



oo 




(3.2) 



n=— oo 




(3.3) 
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But now /a is a characteristic function, so is or 1. It follows that 

exp (-ik6tf A ) = / A + e- iMt /A (3.4) 
where /a = 1 — /a is the characteristic function for the region A outside A, and therefore 

/rib °° 

J] [f A (t n ) + e- ikSt f A (t n )] (3.5) 
n=— oo 

This result has a very appealing form. When the product is expanded out, we get sums 

of products of the characteristic functions /a and f A , so the first term, for example, is the 
continuum limit of 

oo 

n /aw (3.6) 

n=— oo 

This quantity is clearly equal to 1 for a classical history in which the particle is outside A 
at every point along its trajectory and is zero otherwise. The other terms involve similar 
histories including the function /a, so these are histories which enter A for some of the time. 
The integration over k produces a ^-function ensuring that only histories which spend time 
less than e in the region A are included. In particular, as e — > 0, the only term that is left 
is the first term, Eq. ([3.6|) . This is reparametrization- invariant because the expression it was 
derived from is. The important conclusion from this is that we might therefore expect to 
obtain reparametrization-invariant class operators in the quantum theory by taking infinite 
products of projection operators. 

Turning now to the quantum theory, it is well-known that operators commuting with H 
can be constructed using the formula, 



A= dt B(t) (3.7) 

J — oo 

where 

B(t) = e im Be~ iHt (3.8) 
The operator A commutes with H because 

/•OO /»oo 
e iHr Ae -iHr / dt B(t + t) = I dt B(t) = A (3.9) 



so 

[H,A] = (3.10) 
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Suppose we let B = In b. Then, very loosely speaking 

/DC / 00 \ 

dt B(t) = In I [ &(*) (3.11) 
00 \i=-oo / 

That is, to the extent that the continuous product over t is defined, we expect that operators 
of the form 

00 

n & w ( 3 - 12 ) 

will commute with H . 

Given these motivational remarks, we now give the new proposal for class operators for 
trajectories that never enter the region A. As before, denote by P the projector onto A 
and P the projector onto the outside of A. Then our proposal for the class operator for 

trajectories not entering A is the time-ordered infinite product, 

00 

ca= n p{ f) ( 3 - 13 ) 
t=~ 00 

To define this more precisely, we first consider the product of projectors at a discrete set 
of times, i! = t < t\ < £ 2 < • ■■tn-i < t n = t", where tk+i — tfc = St. We define the 
intermediate quantity, C&(t",tf) as the continuum limit of the product of projectors, 

C A (t",t') = \imP(t n )...P(t 1 )P(t ) (3.14) 

ot — >0 

where the limit is n — > 00, <5t — > with t" — t' fixed. Finally, the desired class operator is 

C A = lim C A (t",t') (3.15) 

t"—*oo,t'— >~oo 

This new class operator is clearly closely related to the restricted propagator defined 
above, Eq. (|2.11|) but differs by the presence of unitary evolution operators at either end. In 
particular, we have 

C- A {t",t') = e im " g r {t",t')e- im ' (3.16) 

and therefore 

C A = lim e lHt " g r (t",t') e~ iHt ' (3.17) 

which is the most useful form of the class operator. This is the main result of this section. 

As required, the new class operator commutes with H. This is implied by the construc- 
tion, but more explicitly, we have from Eq. (J3.16|) 

e iHs C A (t", t')e- iHs = e iH(t " +s ^ g r (t", f) e' 1 ^ 1 '^ (3.18) 
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This becomes independent of s as t 



n 



oo, t' — > — oo, hence 



[H,C A ]=0 



(3.19) 



Note that there is no reason at this stage why one should not use a different operator 
ordering of the projectors in Eq. (|3.14j) . (This issue will become significant in the bound case 
treated later). Here, we investigate the consequences of the chosen ordering, which appears 
to be the simplest, but keeping in mind that a different choice may be appropriate. 

The class operator Eq. (J3.17|) is quite different from the original proposal for this class 
operator, Eq. (j2.5|) . in that it does not involve an integral over parameter time. Furthermore, 
unlike Eq. (j2.5|) . the new class operator C^(x",x / ) defined in this way does not in general 
vanish when either end point is in A, so it is not perfectly localized in A. In some sense, it 
corresponds to paths which do not enter the region A but are allowed to enter it at infinite 
parameter time. On the other hand, the new class operator is thoroughly compatible with 
the constraint equation, since it commutes with H, whereas Eq. (|2.5jl does not quite satisfy 
the constraint. 

Generally, in the quantum theory, because the position operator does not commute with 
H, there is an incompatibility between localization in configuration space and the constraint 
equation. It is therefore necessary to make a choice as to which of these two requirements 
should be given precedence. The original proposal Eq. (j2.5|) has exact spatial localization, 
but is not fully compatible with the constraint. The new class operators are fully compatible 
with the constraint but are not perfectly localized in configuration space. Hence the current 
approach gives precedence to the constraint equation over localization. 

Also, as noted earlier, the symmetry of reparametrization invariance is quite subtle in that 
Eq. fl2.5p can be argued to be invariant under the configuration space form of reparametriza- 
tions, even though it is not fully compatible with the constraint. The symmetry generated 
by H is slightly larger than the configuration space form of reparametrizations, so in the 
new class operators we are demanding a slightly more restrictive notion of invariance than 
in Eq. f|2.5jl . It would be of interest to explore these subtle differences in greater detail. 

The class operator Ca for entering the region A is now simply defined by 



A more enlightening formula for it may however be obtained using the path decomposition 



Ca = 1 — C A 



(3.20) 
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expansion, Eq. (|2.13|) . In particular, we clearly have 

C A = lim e iHt " g A (t",t') e- lHt ' (3.21) 

t"— >oo,t'— >— oo 

where g&{t",t') is defined by 

g A (x",t"\x',t') = (x"\g A (t",t')W) (3.22) 

and the left-hand side is given by the path decomposition expansion, Eq. (j2.13|) . 

It is not immediately clear from the definition of these class operators that they will exist 
in all situations of interest. In particular, one would expect that the continuous products over 
time and infinite limits will require careful attention. The proper mathematical framework 
for handling these quantities is the continuous tensor product structure defined by Isham 
et al. Here, we will proceed in a more informal way, and we will see by explicit 

computation in specific examples that the class operator exists and gives reasonable results. 
A more rigoro us ap proach to quantizing models of this type is being pursued by Anastopoulos 



is ap 
Q, 



and Savvidou |13| . using the structures developed by Isham et al. j3J|. Future papers will 



address the connection between the present approach and these more rigorous approaches. 



IV. COMPARISON WITH THE EVOLVING CONSTANTS METHOD 



Willi IIICCVUIVIIJ 



The decoherent histories approach considered here for timeless theories bears comparison 
with the evolving: constants method of Rovelli |l| (for further developments see Refs. 

In that method, one constructs operators commuting with the constraint 
corresponding to physically interesting questions. For these operators, one may construct 
projections P a onto ranges of the spectrum and the probabilities then have the usual form 
Tr(P a p). 

For example, suppose the system is a free particle in two dimensions, with Hamiltonian 

H = A + A ( 4.i) 

2m 2m y 1 

Suppose we are interested in the question, what is the value of x\ when x 2 = t? The 
corresponding evolving constants variable is 

X 1 (r) = ds Xl (s) 5(x 2 (s) - r) (4.2) 

J-oo ds 
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where 



Xi(s) = Xi + —s (4.3) 
m 



for i — 1, 2. This clearly has vanishing Poisson bracket with H. Classically, the integral over 
s may be carried out with the result 

X l (r)=x 1 + ^(r-x 2 ) (4.4) 

P2 

The l/p 2 factor presents difficulties in turning this into a self-adjoint operator, and as a 
consequence the spectrum of states is not orthonormal. One can still construct a POVM 
onto a range of the spectrum but it will not be an exact projector, since it will not satisfy 
PaPp = for a (3. This leads to a kind of imprecision in their definition. These POVMs 
however, may still be useful, in the same way the phase space localized quasi-projectors are 
useful. 

To compare with the decoherent histories approach, suppose we are interested in the prob- 
ability of entering or not entering a region A. Consider therefore the expression Eq. (j2.15|) 
for the parameter time spent in A, which in this simple model is 



r A = ds f A (x 1 (s),x 2 (s)) (4.5) 

J — oo 

Eq. (j4.5p may be written 

t a = ds dy 1 dy 2 5(xi(s) - yi) 5(x 2 {s) - y 2 ) (4.6) 

J-oo J A 

The s integral may then be done with the result 

f Tfl P\ 

ta= dy 1 dy 2 — 5(x 1 ^ (y 2 - x 2 ) - y x ) (4.7) 

J A P2 P2 

Importantly, the result depends only on the evolving constants operator Xi(y 2 ) and on p 2 
(both of which commute with H). Hence in the evolving constraints approach one would 
consider the spectrum of the operator r A , using what is known about the spectrum of Xi(y 2 ) 
and p 2 , and attempt to construct a projector or POVM onto ranges of the spectrum of r A 
(bearing in mind the difficulties noted above with self-adjointness). We can then find the 
probability of not entering A using the quasi-projector 6(e — t a ). 

In the decoherent histories approach, we also take Eq. (|4.5|) as the starting point. However, 
from that we deduce the classical expression Eq. (|3.6j) . which may be used for computation 
of the probability of not entering A. This is the starting point for the quantum theory, 
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and in particular, it inspires the construction of class operators in terms of products of 
projection operators, as described in the previous section. Importantly, class operators are 
not required to be self-adjoint operators, which in some sense means there is more freedom in 
the decoherent histories approach. On the other hand, in the decoherent histories approach 
probabilities cannot be defined in general, unless there is decoherence, so in this sense the 
theory is more restrictive than the evolving constants method. 

The two resulting quantum theories are clearly quite different. However, what they have 
in common is that they take starting points which are classically equivalent. A more detailed 
comparison of these two approaches will be undertaken elsewhere. 

V. THE PARAMETRIZED NON-RELATIVISTIC PARTICLE 

When developing a quantization scheme for parametrized systems, one of the most im- 
portant simple systems to apply it to is the parametrized non-relativistic particle. This is 
because its quantum theory is standard non-relativistic quantum mechanics and it is there- 
fore easy to check whether the expected results are reproduced by the methods described in 
Section 3. 

A. The Parametrized Particle and its Quantization 

The parametrized particle is the usual non-relativistic particle but with the time coordi- 
nate t raised to the status of a dynamical variable, with conjugate momentum p t . Its action 
in Hamiltonian form is 



where a dot denotes differentiation with respect to the unphysical time parameter s. (Note 
that t is physical time in this section). N is a Lagrange multiplier enforcing the constraint 




(5.1) 



H = p t + h = 



(5.2) 



where h is the usual Hamiltonian 



Canonical quantization leads to the Schrodinger equation 




(5.3) 




(5.4) 
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We are ultimately interested in solutions to the constraint equation, Eq. (j5.4j) . which are 
normalized in terms of an inner product defined on spacetime (not just on space). Following 
the general scheme for constructing the induced inner product, we consider an enlarged 
Hilbert space Ti x ®Ti t , where Ti x is the usual Hilbert space of wave functions ip(x). We may 
define states on this enlarged space of the form 

= / dxdt \x) ® \t) V(x,t) (5.5) 



where \x) and \t) are eigenstates of the operators x and t respectively. We then consider 
eigenstates of H , 

H\%) = X\%) (5.6) 
They are normalized using the auxiliary inner product defined on Ti x ® Ti t , 

(#,|#{,) A = J dxdt^{x,t)%(x,t) (5.7) 
Since H = p t + h, the solutions to the eigenvalue equation may be written 

%(x,t) = -I-? e at 1>(x,t) (5.8) 

(27r) 2 

where ip(x,t) satisfies the Schrodinger equation, Eq. (j5.4|) . It follows that 

(*l\%) A = ^j dt j dx e ~ U+lVt r{x,t)iP'{x,t) (5.9) 

The integral contains within it the usual inner product 

= J dx ip*{x,t)ij)'{x,t) (5.10) 

This has the important property that it is independent of time when the states obey the 
Schrodinger equation, so the time integral may be done in Eq. (j5.9|) . pulling down a delta 
function 6(1 — I'). We thus obtain 

(%\%) A = 6(l-l'){iPW) s (5.11) 

This means that the expected Schrodinger inner product on surfaces of constant t is fully 
compatible with the induced inner product defined on the whole of spacetime. 

We may now construct the decoherence functional for this system, for some interesting 
physical questions. For a pure initial state, the decoherence functional is 

D(a,a') = (%,\ClC a \%) (5.12) 

It is constructed using the induced inner product and the class operators C a must commute 
with the constraint, H. 
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B. A Useful Mathematical Result 



The induced inner produce we are using here can be quite cumbersome and it is useful 
to prove a simple mathematical result. The decoherence functional is an expression of the 
general form 

(VvlAlVt) (5.13) 

where A commutes with the constraint, and it is useful to show how this expression reduces 
to a simpler expression on the original Hilbert space Ti x . In all the expressions we are 
interested in, A will commute with t, so A has the form 

/oo 
dt B(t) ® (5.14) 
-oo 

where B(t) acts on Tt x only. Using the fact that [A, H] = 0, it is straightforward to deduce 
that 

B{t) = e- iht B{0)e lht (5.15) 

(where note that the signs in the exponents are not the ones associated with Heisenberg 
picture evolution). We now have 



= J dt J dx'dx %(x',t){x'\B(t)\x) %(x,t) 



2vr 

Now noting that 



— / dt e - ilt+m I dx'dx i)*{x',t) (x'\e- iht B(0)e m \x) ij){x,t) (5.16) 



dx e iht \x) *l){x,t) = I dx \x) ip(x,0) = |V>) (5.17) 



we finally obtain 

(* V \A\%) = 6{1-1') (V|S(0)|V>) (5.18) 

The expression on the left is in terms of the auxiliary inner product on Ti x <g> Tit- The inner 
product on the right is the usual one on Ti x . Hence expressions of the form Eq. (|5.13J) are 
readily evaluated once one has read off B(t) in Eq. (j5.15j) . (Note these expressions specifically 
refer to the parametrized non-relativistic particle with Hamiltonian constraint Eq. (|5.2j) - 
they are not valid for parametrized systems which are quadratic in all the momenta). 
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C. Probabilities on Surfaces of Constant Time 



Now we consider the simple question, what is the probability of finding the particle in a 
range A of the x-axis at time t^l We of course expect the standard answer 

PA= I dx \ip(x,t )\ 2 (5.19) 

J A 

but it is important to see how this arises in a consistent quantization of the parametrized 
particle. 

We assert that the appropriate class operator corresponding to this question is 



C A ,t = / ds 6(t(s) - t ) f A (x(s)) (5.20) 

J — oo 

where 

t(s) = e iHs te- iHs = t + s (5.21) 

x (s) = e iHs xe~ iHs (5.22) 

and Ja{x) is a window function on the range A. The class operator clearly commutes 
with both H and t. Classically, for the free particle, the class operator corresponds to the 
expression 

C A ,t„ = /A(a:--(t-to)) (5.23) 
m 

which clearly has the right properties: it is equal to 1 for classical trajectories which cross 
t = to in A and zero otherwise. Returning to the quantum case, one can see that 
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C A ,t = / dt e- ih(t - to) f A (x)e ih{t - to) ®|t)(t| (5.24) 

This is of the form Eq. (|5.14|) from which we read off 

B(0) = e iht0 f A {x)e- iht0 (5.25) 

Now a crucial simplification. Since /a is a window function, B(0) is in fact a projection 

operator, and therefore the class operator C A ,t is also a projection operator. This means 
that decoherence is automatic, between histories characterized by C^ to and 1 — C& )to and 
we may immediately assign the probability 

PA = <#HCA,tol*i) (5-26) 
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Using Eq. (|5.18p and dropping the 5- function, this becomes 

PA = (iP\e m °f A (x)e- iht0 \<iP) (5.27) 

which agrees exactly with the expected result Eq. (|5.19Jl . Note also that since the class 
operator is a projection operator in this case, the decoherent histories analysis agrees exactly 
with the evolving constants approach. 



D. Probabilities for Spacetime Regions 

We now consider a more challenging question which is to consider probabilities for regions 
of spacetime. In particular, we pose the following question: given the initial wave function 
ip{x) at t — 0, what is the probability that the particle is found in the region x < in the 
time interval [0, r]? This has been analyzed previously in the decoherent histories approach 
with the following results (See also Refs.j^J, S])- The decoherence functional is 

D{a,af) = (if>\glg a \i[>) (5.28) 

in the usual inner product. (Here we use g a to denote class operators to avoid confusion 
with the class operators C a defined on the enlarged space.) There are two class operators. 
First, there is the class operator for remaining in x > during the time interval [0, r], and 
this is the restricted propagator g r (r, 0). The other class operator may be written 

g c (T,0)=e- ihT -g r (T,0) (5.29) 

and corresponds to the situation in which the particle enters x < at some time during [0, r]. 
(It may also be expressed in terms of the path decomposition expansion Eq. (j2.1Hj) . but we will 
not need this here. Note also that these class operators reduce to the unitary operator e~ %hT 
when the restrictions are removed). The histories are generally not decoherent, except for 
very special initial states, and the resultant probabilities are somewhat trivial 33| . However, 
our aim here is to show how the decoherence functional for this model is recovered from a 
quantization of the parameterized particle as a constrained system on an enlarged Hilbert 
space, in which its spacetime character is most transparent. 

We take the decoherence functional Eq. (j5.12|) and seek a class operator commuting with 
H = p t + h corresponding to the statement that the particle never enters the region x < 
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during the time interval [0, r]. We denote this region A and we use A to denote the region 
outside A. In the geuninely spacetime point of view used here, we may introduce projection 
operators onto the spacetime regions A and A. The projection onto A is 

P = 6( T - t)6{t)6(-x) (5.30) 

and the projection onto A is conveniently written 

P = 6(-t) + 6{t - i)6(t)6(x) + 6{t - r) (5.31) 

/oo 

dt T(t,x) ® \t)(t\ (5.32) 
-oo 

where T(t, x) is an operator on TC X , equal to 8(x) for < t < r and equal to the identity 
otherwise. The class operator for remaining in A (that is, never entering the region A) is of 
the form 

oo 

C r = ] [ P(s) (5.33) 



s=oo 



This example is sufficiently simple that we can work directly with the infinite product (time- 
ordered) without encountering difficulties. We have, from Eq. (|5.32|) . 

/oo °° 
dt ] [ T(t + s,x(s)) ® \t)(t\ 
■°° s=oo 

/oo ZZll 
dt l\ 9(x(s))®\t)(t\ (5.34) 
■°° s=-t 

From the definition of the restricted propagator, Eq. (|2.11Jl . we see that 

/oo 

dt e ih{ - T -Vg r {r, 0)e iht ® |*> <*| (5.35) 
-oo 

This commutes with the constraint H and is of the desired from Eq. (|5.14jl with 

B(0)=e ihT g r (r,0) (5.36) 
We are also interested in the quantity 

/oo 
dt e- iht g r (r,0yg r (T,0)e iht ® \t)(t\ (5.37) 
-oo 

which is also of the form Eq. (j5.14j) with 

B(0)=^(r,0)V(r,0) (5.38) 
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From these objects one can also construct the class operator for crossing, 



(5.39) 



and related objects such as C\C r . It is now easy to see, using Eq. (j5.18j) . that we readily 
obtain the known form of the decoherence functional for this system. For example, 



which, via the induced inner product prescription, agrees with the known result Eq. ([5.28)1 . 

These results show that our proposal for class operators passes the important test of 
the quantization of the non-relativistic particle in parametrized form. Furthermore, there 
is the added feature that it shows how spacetime questions in non-relativistic quantum 
mechanics may be expressed in a genuinely spacetime form, since the decoherence functional 
and probabilities (such as the left-hand side of Eq. ([5.40)1 ) may be expressed in terms of an 
inner product and operators defined on spacetime. 

VI. A SIMPLE ONE-DIMENSIONAL EXAMPLE 

The parametrized non-relativistic particle has the very special simplifying feature that 
the Hamiltonian is linear in one of the momenta. This is not the case in general. We 
therefore now consider some examples with a Hamiltonian quadratic in the momenta. We 
first consider a simple one-dimensional example involving the free particle. It is trivial in 
itself (except to show that the class operators can be easily calculated), but readily extends 
to higher dimension and has important implications for the relativistic particle considered 
later. 

A. Energy Eigenstates 

We first consider normalization of the energy eigenstates. We have 



where H = p 2 /2m. For each E there are two solutions which are conveniently written, 




(5.40) 




(6.1) 




(6.2) 
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where \k\ = \/2mE. In the usual inner product, 

(H^2) = J dx il>\{x)il) 2 (x) (6.3) 

we have 

{^%) = 8{E - E') (6.4) 

and 

<^il^> = (6.5) 

In the induced inner product prescription we therefore drop the 5-function and take the 
induced inner product between two eigenstates with the same E to be 

<^I>7 = 1 (6-6) 

and 

(^|>7 = (6.7) 



B. Class Operators for Crossing or Not Crossing the Origin 

Given these preliminaries, we now consider the following simple question. Given that the 
system is in an energy eigenstate, what is the probability that the particle crosses or never 
crosses x — 0, irrespective of time? This is most easily handled by considering the class 
operator C r describing the situation in which the particle is always in x > or x < 0. The 
class operator C c for crossing x = is then given by 

C c = 1 - C r (6.8) 

Let P be the projector onto the positive x-axis, 

poo 

P= dx \x)(x\ = 6{x) (6.9) 
Jo 

(we use the hat symbol for operators only where clarity demands it). Then 

P=l-P = 6(-x) (6.10) 

is the projector onto the negative a;-axis. The class operator C r for remaining in x > or 
x < is then given by, in a loose notation, 

Cr = II P W + II ^(*) = C r + C; (6.11) 

t t 
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where + and — denote the terms projecting onto the positive and negative x-axis respectively. 
This expression is defined more formally in terms of restricted propagators, as in Eq. (j3.17|) . 

This situation is sufficiently simple for the method of images to work. The restricted 
propagator for the region x > is therefore 

g+(x", t\x', 0) = 9{x")9{x') [g(x", t\x', 0) - g(x", t\ - x', 0)] (6.12) 

where 

g(x",t\x',0)= (^^expfe^-^ 2 ) ( 6 - 13 ) 

is the free particle propagator. The restricted propagator may be usefully thought of as sum 
of two parts. The first term (in a semiclassical view) corresponds to the direct path from 
x' to x" . The second term is usually thought of as propagation from the image point —x'. 
However, it may also be thought of as corresponding to the path which again starts at x' 
but is reflected off x = before arriving at x" . Differently put, the classical limit of a system 
described by a restricted propagator is one in which the Hamiltonian includes an infinite 
potential barrier at the boundary of the region in question, so its classical trajectories include 
paths which reflect off the boundary. These points will be important in the interpretation 
of the quantum results. 

It is in fact quite useful to rewrite this in an operator form using the projection operators 
introduced above. We also introduce the reflection operator 

/oo /*oo 

dx \x)(—x\ = dp \p)(—p\ (6.14) 
-oo J — oo 

and note that [H,R] = and R 2 = 1. The restricted propagator in x > may then be 
written in operator form as 

g +(t, 0) = Pe~ iHt (l - R)P (6.15) 
From Eq. (|3.17jl . the desired class operator is 

C+= lim P(t") (1 - R) P(t') (6.16) 

i"^oo,t'^— oo 

Now we need to take the infinite time limit in P(t") and P(t'). We have 

P(t)=9(x t ) = e(x + ^) (6.17) 

m 
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Naively, for very large positive or negative t the momentum term dominates, and we get 



lim P{t") = 6(p) (6.18) 

t"^oo 

lim P(t') = 0(-p) (6.19) 

f— >— oo 

However, it is important to verify this result more carefully, by sandwiching it between an 
arbitrary pair of states. We have, using Eq. (j6.13J) . 



m r°° poo poo 
(il)i\P(t)\i/) 2 ) = — J dx J dy 1 J dy 2 ipKyx) il) 2 (y 2 ) 



J —oo J —oo 

im , 2 im 



x exp ( -—{ Vl - xY + —(x - y 2 y ) (6.20) 



Defining p = mx/t, this becomes 



<Vi|P(t)|V2> = ^ / <IP I d m I d, l2 r :(//,) r 2 (// 2 ) 



-oo 

im 



x exp lip( yi - y 2 ) - ^-{yl - y 2 2 ) j (6.21) 

Now, taking the limit t —>■ oo, the last term in the exponent drops out, and the integrals 
over y\ and y 2 produce the Fourier transformed wave functions, so we have 

roo 

lim(^i|P(t)|^ 2 ) = / dp ipl(p) ip 2 (p) 

t^-oo 







(HmM (6.22) 



This confirms the naive result Eq. ()6.18J) . 

We now have 



which may also be written 



C+ = 9{p){l - R)9{-p) 

= -6(p) R 0(-p) (6.23) 

C+ = -9{p)R = -R6(-p) (6.24) 



Note that the first part of the restricted propagator (involving the direct path from initial 
to final point) drops out, leaving only the part corresponding to the reflected paths. In the 
momentum representation, this is 

(k 2 \C+\h) = -6(k 2 )5(k 2 + h) (6.25) 
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Similarly, the class operator C r is constructed using the restricted propagator in x < 0, 
g~(x", t\x', 0) = 6(-x")6(-x') [g(x", t\x\ 0) - g(x", t\ - x', 0)] (6.26) 
and one readily finds that 

C; = -6(-p)R = -R6(p) (6.27) 
Combining the + and — parts we have 

C r = -G(p)R - 0(-p)R = -R (6.28) 

a very simple expression. The class operator for crossing the surface is then 

C c = 1 - C r = 1 + R (6.29) 

As expected, all the above class operators, C c , C r , and C~ commute with H. 

C. Decoherence Functional and Probabilities 

We may now compute the decoherence functional and the probabilities. The off-diagonal 
part of the decoherence functional is 

D(r,c) = Tr (C r pCi) (6.30) 

computed in the induced inner product. We take a pure initial state which is an eigenstate 
of H. Noting that 

ClC r = -(l + R)R = -(R + l) (6.31) 

we have 

D(r,c) = -(* E \(1 + R)\* E ) I (6.32) 

where \^e) is an energy eigenstate and the induced inner product is used. Using the formula 
Eq. (jl.21j) for the induced inner product, this means that there is decoherence only for states 
satisfying 

((j)\5{H - E){l + R)\<j)) = (6.33) 

where \^e) — 8{H — E)\<p). (This expression is in fact real, so there is no difference 
between consistency and decoherence). The operator (1 + R) produces the symmetric part 
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of The decoherence condition therefore implies that the fiducial wave function must be 
antisymmetric, 

(j)(-x) = -<f>(x) (6.34) 

This means that ^e{x) must also be antisymmetric (using the fact that R commutes with 
H). For such wave functions, probabilities are defined and we get 

Pr = ^ E \ClC r \m E ) I = l (6.35) 

for the probability for not crossing since C\C T — R 2 — 1. Similarly, the probability for 
crossing is zero. 

It is useful to try and understand this result in terms of classical paths. Recall that 
firstly, only entire infinite classical paths are reparmetrization invariant and secondly, that 
the restricted propagator corresponds to a classical situation in which there is an infinite 
barrier at x — 0. Classically, therefore, this result in some sense corresponds to classical 
trajectories which remain in x > or x < by bouncing off x = 0. 

This result is not very physically enlightening but it is very similar to the result obtained 

in the decoherent histories analysis of the crossing time problem in non-relativistic quantum 

mechanics. There, one looks for probabilities that, given an initial state, the particle will 

cross x = during the time interval [0, r\. One finds that only antisymmetric wave functions 

give consistency and the crossing probability is zero [3] . More physically intuitive results in 

the crossing time problem are obtained — when there is a decoherence mechanism in place 
~l 

|. We expect that to be the case here, too, but this will explored in another paper 



D. A More Detailed Look at the Crossing Class Operator 

It is useful to give an alternative derivation of the crossing class operator Eq. ()6.29|) using 
the path decomposition expansion, Eq. (j2.13J) . This is partly a consistency check but it will 
also give some insight into the form of the result. 

We first consider the fixed time crossing propagator, taking into account the possibility 
of crossing the origin in either direction. Applying Eq. (J2.13|) together with the restricted 
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propagators Eqs. (|6.12|h ()6.26|) . this is 



r 1 " % 8 

g c (x",t"\x',t f ) = lim / dt g(x", t"\x, t) — — g+(x, t\x', t') (6.36) 

f % d 

— lim / dt g(x" ,t"\x,t) - — —g r (x,t\x',t') 



x^o+ J t , ' ' 2m dx * 

where the relative minus sign is because of the definition of the normal in Eq. (|2.13jl . Now, 
using Eqs. (j6.12j) . ()6.26j) . we have 

Jim ^g;(x } t\x',f) = 29{-x') ^-g( x ,t\x',t')\ x=Q (6.37) 

^Mx',t') = 26{x') ^9(x,t\x',t')\ x=0 (6.38) 

We therefore have 

g c (x",t"\x',t') = lim / dt g(x", t"\x, t) — — g(x, t\x', t') e(x') (6.39) 

^ — >0 J y Tfl OX 

where e(x) is the signum function. This is conveniently written in operator form as 

g( t \ t') = - / dt e -^C*"-*) Six) v e~ iH ^ e(-x) (6.40) 
m J t i 

The desired crossing class operator is now given by 

C c = lim e +iHt "g c (t", f) e~ im 

1 
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dt5(x t ) \p\ (6.41) 

where we have used the fact that e(—x t ) — > e(p) as t — > — oo and \p\ = p e(p). 

Eq. (j6.4ip is the desired expression and shows very clearly that the class operator for 
crossing x = involves some kind of flux at x — 0. Classically, it is easy to see that this 
expression is equal to 1 for all classical paths, except those for which p = 0, in which case 
it is zero. Hence it clearly encapsulates the classical notion of surface crossing. One might 
argue that classical states with p = are a set of measure zero so may be safely neglected. 
However, the p = states seem to be crucial to understand the quantum case. 

As an operator expression, Eq. (j6.41|) may be evaluated by sandwiching it between two 
momentum states: 

(p'VcW) = ^J°° dt 6X p(|^ ,/2 -^ 2 )) WW} \p'\ 

= 2\p'\ 5(p" 2 -p' 2 ) 

= 5{p" -p') + 5{p" + p') (6.42) 
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It follows that 



C c = 1 + R 



(6.43) 



the expected result. 

How are we to understand this result? Classically, we noted that all trajectories with p ^ 
cross the origin. The result Eq. (jfi.43|) ought therefore to be the quantum implementation 
of this idea. The key thing is that the operator 1 + R is zero when acting on states which 
are antisymmetric in x about the origin. Such states are also antisymmetric in p so clearly 
vanish at p — 0. Hence the crossing class operator annihilates a class of states with p = 
and in this sense implements the classical notion of crossing. 

Of course, there are many inequivalent ways to turn classical expressions into quantum 
operators, and one could imagine that a quantization procedure may exist which consistently 
drops the p = states before quantization. This would avoid the difficulties of interpretation 
with reflected paths encountered earlier. However, this does not appear to be the case in 
the present quantization method. 

We also remark that, as noted earlier, the evolving constants method also encounters 
difficulties with p = states because of the 1/p factors arising in the evolving constants 
operators. 

VII. THE RELATIVISTIC PARTICLE 

The result of the previous section is readily extended to the relativistic particle in 3 + 1 
dimensions described by the constraint equation 



We consider this in order to compare with a previous attempt to define class operators 
compatible with the constraint |H fl(3|. 

Consider the following question. What is the probability that a free relativistic parti- 
cle never crosses the spacelike surface x° = 0? Classically, this probability must be zero, 
because every (timelike) classical trajectory crosses any spacelike surface somewhere, since 
the classical trajectories are just straight lines. In Ref.[l0(, a (somewhat ad hoc) proposal 
was made to define class operators compatible with the constraint which coincide with this 




(7.1) 
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classical intuition. A quantum-mechanical probability of zero was thus obtained for not 
crossing the surface. 

The question is readily addressed in the present approach using an elementary extension 
of the results of the one-dimensional model in the previous section. The class operator for 
not crossing x° = is of the form Eq. (JbMlj) where the Hamiltonian is as in Eq. (|7.1|) and the 
projectors P and P project onto the regions x° > and x° < respectively. One easily finds 
that decoherence is only possible for states antisymmetric about x° = and the probability 
for not crossing is equal to 1. It is therefore the exactly opposite result to that obtained in 
Ref.Q 

What is the origin of the difference in results and which is the "correct" one? The key 
point is that in the present approach, the restricted propagators involved in the construction 
of the class operators, such as Eq. (j6.12j) involve two types of paths in a semiclassical picture, 
the direct paths and the reflected paths. As noted above (after Eq. (|6.24Jl ). the part of the 
propagators corresponding to the direct paths drops out, so the present approach consists 
entirely of the contribution from the reflected paths. The classical intuition of Ref.[10] was 
tacitly based on the direct paths, so the result is completely different. Since the reflected 
paths capture the important notion of the reflection of wave packets from a barrier, it is 
appropriate to take the present approach as the definitive one if we are to stay true to 
quantum-mechanical ideas. 

In the closely related context of the arrival time problem, it has however been shown that, 
in the presence of a decohering environment, the effect of the reflected paths becomes ess 

n n 

significant, and this is how classical intuition may become restored |3JJ. See also Ref. |37| 
for further relevant considerations of the relativistic particle. 

This brief discussion of the relativistic particle indicates that the present proposal for class 
operators is not in fact a more developed statement of the ad hoc approach of Refs. 0, [it], [ll | , 
but a different proposal altogether. 



VIII. THE FREE PARTICLE IN TWO DIMENSIONS 

In two dimensions more interesting questions are possible. However, the general difficulty 
one expects to encounter for most questions is decoherence. In the absence of an environ- 
ment, most situations will not have decoherent sets of histories. In the one-dimensional 
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example, the lack of decoherence for general states is largely due to the feature of reflection 
at the boundaries of the regions of interest. It is therefore of interest to consider situations 
where this reflection will not happen. 

In two (and more) dimensions, we have the possibility of eigenstates of H which are 
rotationally symmetric about the origin and may be thought of as superpositions of wave 
packets moving radially. It therefore seems plausible that if we consider questions concerning 
regions whose boundaries lie along radial lines, there will be little or no possibility of crossing 
or reflection (since the wave function has no flux across the boundary). For example, given 
a rotationally symmetric state, we could ask for the probability that the particle is found 
in a wedge-shaped region emanating from the origin. This question bears some resemblance 
to Mott's calculation of alpha-ray tracks He showed, using a series of model detectors, 
why an outgoing spherical wave produces a straight line track in a detector. 

A. Particle in The First Quadrant 

We first consider the following question. Given that the system is in an energy eigenstate, 
what is the probability that the particle will always be in the region, x > 0, y > 0? We 
denote this region A and use the method of images to construct the restricted propagator 
for propagation in A. It is 



where g(x", y", t\x', y', 0) is the free particle propagator in two dimensions. As in the previous 
example, an operator form is useful. We introduce two reflection operators 



g A (x"y,t\x',y',0) = d{x")e{y") 6(x f )6(y f ) 

x [g(x", y", t\x', y', 0) - g(x", y" , t\ - x', y\ 0) 
- g{x", y", t\x', -y', 0) + g{x", y", t\ - x', -y', 0)] 



(8.1) 




(8.2) 




The class operator for histories which are always in this region is therefore 





(8.4) 



oo 
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where 



P(t) = 0(x t )9(y t ) 



(8.5) 



Following the same method as in the previous example, we therefore get 



C A = e{p x )6(p y ) (1-R x )(l-R y ) 6{-p x )6{-p y ) 



(8.6) 



Using the properties of the reflection operators this is easily seen to be, 



C A = 9{p x )9{p y )R x R, 



■y 



(8.7) 



For a pure initial state, the off-diagonal part of the decoherence functional is 



D(A, A) = (m-C A ^C A \iP)j 

= tyKRtRy-lM-pjei-PvMi 



(8.8) 



where \ip) are energy eigenstates and the induced inner product is used. (In what follows 
we will not spell out the induced inner product calculation since it is very similar to that in 
Section 6). From this it is easy to see that we get decoherence if the state satisfies 



This condition is easily satisfied, as expected, by a rotationally symmetric state and the 
probability for remaining in the region A is 



For rotationally symmetric states this will be equal to 1/4, the expected result. 

Classically, this result corresponds to straight line trajectories radiating from the origin 
in the positive quadrant. One can think of them as coming in from infinity, bouncing at the 
origin and returning to infinity. 

B. Particle in the First and Third Quadrant 

Since the classical trajectories are straight lines, as a modification of the above calculation, 
it seems reasonable to also consider the situation in which the region A is the first and third 
quadrant. We can then ask whether it is possible to have a situation which corresponds 
classically to a family of infinite straight lines passing through the origin. 





PA=(4>\9(PMPyM)l 



(8.10) 
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There is a subtlety, however, in that there are two similar but different class operators. 
First, one can ask for the probability that the particle remains always in the first quadrant 
x > 0, y > or always in the third quadrant, x < 0,y < 0. It is easy to see that the 
appropriate class operator has the loose form 

c A =n p iw+n p 3w (8.H) 

where Pi and P3 are projectors onto the first and third quadrants, respectively. Or more 
formally, this is 

C A = [0(p x )0(Py) + 0(-p x )0(-Py)] RxRy (8.12) 

We again get decoherence for states satisfying Eq. ()8.9|) and the probability associated with 
A is 1/2. 

Second, one can ask for the probability that the particle is in the first or third quadrant, 
but in addition, has the possibility of following trajectories which lie in both quadrants. 
(This possibility is not present in the class operator Eq. (j8.12j) ). The class operator for this 
second case has the loose form, 

C A = H(P 1 (t) + P 3 (t)) (8.13) 

t 

which is clearly different from Eq. (j8.11|) . By comparing with Eq. (j8.4|) . it is easily see that 
this class operator is 

C A = [OiPMPy) + 0(-pM-Pv)) (! " R * ~ R v + R * R y) WxWv) + 0(-P*)0(-Pv)] 
= [6(p x )6(Py) + d(-PM-Py)] (1 + RxRy) (8-14) 

It is easy to show that 

c1c a = 2Ca (8.15) 
from which it follows that the off-diagonal term in the decoherence functional is 

D(A,A) = (^KI-Ca^CaI^j 

= -MC a |V)j (8.16) 

The factor of (1 + R x R y ) in C A means that there is only decoherence for states satisfying 

R X Ry\lp) = ~\lp) (8.17) 
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This is the opposite condition to the previous case and is not satisfied by the rotationally 
symmetric wave function. Moreover, the off-diagonal part of the decoherence functional is 
in fact proportional to the probability, so for states satisfying the decoherence condition, the 
probability associated with A is zero. There does not appear to be an obvious interpretation 
of this result in terms of classical trajectories. 

C. Particle in a Wedge-Shaped Region 

We now consider the case of a more general region A consisting of the wedge lying in 
the region < <fi < (3 (in polar coordinates r, 0). We again ask for the probability that the 
particle is always in the region A. 

Following the general scheme, we first require the time-dependent propagator for the 
wedge region. For simplicity we restrict to the case where the angle (3 is (3 = 7r/6, where b is 
an integer. We also take b to be even (which turns out to be simplest to deal with). Then 
the restricted propagator for the interior of the region is 



where g(r, <p, t\ro, <po, 0) is the free particle propagator in two dimensions in polar coordinates, 
and fp(x,y) is a characteristic function equal to 1 inside the wedge region and zero outside 
. The desired class operator is now 



gp(x,y,t\x o ,y o ,0) = fp{x,y) fp(x ,y ) 



6-1 



[9(r, 2nf3 + 0, t\r , <f> , 0) - g(r, 2n(3 - 0, t\r , O , 0)] (8.18) 



n=0 





n=0 




where the projector P is f/3(x,y) and we have introduced the rotation operators 





(8.21) 



Using the same method as in Eq. (j6.22J) . it may be shown that 
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so we obtain 

6-1 

C(3 = fp(Px,Py) \- Rn ~ K ^ fp(-P*> -Pv) ( 8 - 23 ) 

n=0 

Since R n rotates by an angle 2n/3, we have 

fp(Px,Py)Rnfp(-Px, -Py) = fp{Px,Py) M~P^ -Py)^ (8-24) 

where p™, p™ denote the momenta rotated through angle 2n(3. Clearly, Eq. (j8.24|) is zero 
unless n = b/2 (recall that b is an even integer). Similarly, it is readily shown that 

fp(Px,Py)K n f /3 (-p x , -p y ) = (8.25) 

We now have 

Cf3 = fp{p x ,P y )R b/ 2 (8.26) 

This agrees with Eq. (j8.7J) when (3 = it/ 2 as expected. 
The off-diagonal term in the decoherence functional is 

= (i/)\{Rb/2 - l)f{-Px, -PvWi (8-27) 

in the induced inner product, where |?/>) are energy eigenstates. This means that there is 
decoherence for states satisfying 

R b /M = (8.28) 
This will indeed be satisfied for the rotationally symmetric state and the probability then is 

pp = MCpMi 

= Wp(Px,PyM)i (8.29) 

By symmetry, we clearly have 

Pfi = ^ (8-30) 

as expected. 

To summarize this section, we have obtained some reasonable physical results for rotation- 
ally symmetric wave functions and wedge-shaped regions, although some simple extensions 
of these ideas run into some interpretational difficulties. These examples uphold some of 
our intuitions, and it is gratifying that it is not always necessary to employ an environment 
to get decoherence. 
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IX. SYSTEMS OF HARMONIC OSCILLATORS 



The formalism so far concerned unbound systems, in which the (unphysical) time- 
parameter t runs from — oo to +00. It is however very different (and simpler) for systems 
of harmonic oscillators, which are periodic in time. In this section we consider the case of a 
(i-dimensional simple harmonic oscillator with Hamiltonian 

H = i(p 2 + x 2 ) (9.1) 

(See Ref.jjJ for the evolving constants analysis of this system). Much of the formalism 
will, however, be applicable to other systems periodic in time. Systems described by the 
Hamiltonian Eq. (j9.1j) will have period 2it so in the quantum theory we have 

el H(t+2ir) = e iHt ( 92 ) 

An important class of observables for this system are of the form 

2tt 



A= dt B(t) (9.3) 
Jo 

It is easy to show that A commutes with H using the fact that the periodicity implies the 
property 

2-7T / t 27T+T 

dt B(t) = dt B(t) (9.4) 



Jt 

for any r. Note that because the spectrum of H is discrete, it is not necessary to use the 
induced inner product. 

A. Class Operators 

For the systems described by Hamiltonian Eq. (|9.1|) . the natural modification of Eq. ()3.13jh 
the class operator for not entering the region A, is 

<7 A = Ip(t) ( 9 - 5 ) 
t=o 

where P is the projector onto the region outside A. However, it is not hard to see that this 
does not in fact commute with H (unless the P(t) all commute at different times). This is 
an operator ordering issue and is easily remedied by defining the class operator to be instead 



•I p2n s + 27r 

Ca = ^ J o ds]} P(t) (9.6) 
37 



G K = i- / e lHs g r (2n, 0) e" lHs (9.8) 
o 



C A = ^- I ds e iHs g c (2rt, 0) e~ iHs (9.9) 



o 



This is essentially a sum over all cyclic permutations of the operators P(t) at different times, 
and now commutes with the Hamiltonian. 

Following steps similar to those used in previous Sections, it is easily shown that 

i r 2n 

C A = — ds e ^ +2w) g r {s + 2tt, s) e~ iHs (9.7) 
27r Jo 

where g r is the restricted propagator. (Note that the restricted propagator will not in general 
be periodic in time.) Using Eq. (j9.2j) together with the fact that g r (t,t') depends on time 
only through (t — t') , we have 

/■27T 

This is of the form Eq. (j9.3|) so commutes with H, as expected. 
The class operator Ca = 1 — C*a rnay be written, 

/■27T 

2tt 

where gd^Tr, 0) is the crossing propagator, given by a sum over paths which enter A at some 
time during the interval [0, 2tt]. 

B. Decoherence Functional and Probabilities 

We may now look at the decoherence functional and the probabilities. We choose a pure 
initial state \ip) and, in keeping with the general approach, this state is taken to be an 
eigenstate of the Hamiltonian 

H\tfj) = E\tfj) (9.10) 
The off-diagonal term of the decoherence functional is 

D(A,A) = (V;|C a Ca|V>> (9.H) 
Inserting the explicit expressions, Eqs. (j9.8|) . (j9.9j) . we have 

IT I ^l^ Sl ^(27r,0)e-^-^^(27r,0)e- iHs ^) (9.12) 

Using Eq. (j9.1U|) this becomes 

D(A, A) = (Vb2(27r,0) P E ^(2tt,0)|V) (9.13) 
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where we have introduced the object 

Pe= r ?± e -i{H- E )s (914) 

Jo 2n 

Because the spectrum of H is discrete, this is a projection operator, so satisfies P| = Pe- 
The decoherence functional will not be diagonal in general, although one simple case in 

which it will is when the wave function is an eigenstate of g c (27r,0) or g r (27T, 0). We will 

exhibit such states below. 

When the decoherence condition is satisfied, the probabilities associated with A and A 

are easily shown to be 

p A = (Vb c (27T,0)|V) 

Pa = (Vk(27r,0)|V> (9.15) 



C. Some Special States Exhibiting Approximate Decoherence 

We now introduce some states which exhibit approximate decoherence and quasiclassi- 
cal behaviour for this model. Consider first the standard coherent states of the harmonic 
oscillator, |p,x). They are preserved in form under unitary evolution, 

e-^|p,x) = e-' i / 2 | Pt ,x t ) (9.16) 

where pt,x 4 are the classical solutions matching p,x at t = 0, hence they are strongly 
peaked about the classical path. There is a set of states which are natural analo gues of 
these states for the timeless models considered here. They were referred to in Ref.|40( as 
"timeless coherent states" and are defined by 



) = Pb|p,x) 

e -i(H-E)t | p;X) 



^ dt 







2tt 

e^'-a* |p t ,x t > (9.17) 



2tt 



They are clearly eigenstates of H with eigenvalue E and are concentrated around the entire 
classical path with initial data p, x. They are not normalized to 1 exactly, but if the initial 
data satisfies E = |(p 2 + x 2 ) then the coherent states |p,x) are approximate eigenstates 



39 



of Pe and the timeless coherent states are then approximately normalized to 1. Further 
properties of these states are described in Ref.J^. 

Now consider a timeless coherent state whose trajectory Pf,x t lies entirely within the 
region A. This region could, for example, be a large rectangular region in configuration 
space. Or it could be a tube following the classical trajectory but broadened out beyond 
the scale of quantum fluctuations. In both these cases, if P is the projector onto the region 
A, then we clearly have 

^|0px) « |0px) (9.18) 

and 

P|0 px ) « (9.19) 

We assert that with this choice of A, the state |0 px ) will give approximate decoherence. 
There are two ways to see this. 

First, from the (informal) expression Eq. ()9.6|) . the result Eq. (j9.18J) together with the fact 
that the state is an eigenstate of H imply that it is also an approximate eigenstate of P(t), 
so will be an approximate eigenstates of the class operator. This means there is approximate 
decoherence. 

Second, and perhaps a little more rigorously, we use the expression for the decoherence 
functional Eq. ()9.13|) . The important thing is to consider the action of the restricted prop- 
agator g r (27r,0) on the state |0 px ), which, via Eq. (j9.17|) . boils down to its action on the 
coherent state |p, x). Restricted propagators are very difficult to calculate for arbitrary 
regions, but their path integral form gives an intuitive picture of their properties. It is 



g r (x", 27r|x', 0) = J Pxexp (i j dt 



27T 

(9.20) 



1-2 1 2 

-x X 



This is a sum over paths x(t) which remain always in the region A and satisfying the end- 
point conditions x(0) = x', x(27r) = x". Suppose that the trajectory pt,x 4 of the coherent 
state |p, x) remains entirely within the region A (and does not approach the boundary). 
Then, when this initial state is attached to the restricted propagator, the path integral will 
be dominated by the classical path with initial data p, x. The path integral will therefore 
be approximately the same as the unrestricted path integral, which means that 

p (2tt,O)|p,x) « |p,x) (9.21) 
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It follows that there will be approximate decoherence and the probability for finding the 
particle in the region A is approximately 1. 

So for these specially chosen regions that entirely contain the trajectory of the timeless 
coherent state we get approximate decoherence and the expected probabilities. Note that 
these heuristic arguments only work for periodic systems in which certain states remain 
coherent. For the systems considered in earlier sections involving an infinite range of time, 
the spreading of wave packets would render such heuristic arguments invalid. 

For most other choices of A, however, there is no decoherence and probabilities cannot 
be assigned without a decoherence mechanism. This will be pursued elsewhere j^. 



X. SUMMARY AND CONCLUSIONS 



We have discussed the issues involved in defining class operators for the decoherent his- 
tories analysis of reparametrization invariant systems and made a specific proposal for such 
operators. The class operators defined are based on certain reasonable classical expressions 
and reduce to projection operators when everything commutes. They commute with the 
Hamiltonian so fully respect reparametrization invariance. They do not, however, exhibit 
the localization properties of their classical counterparts. This is because there is an incom- 
patibility between localization and the constraint and in our definition we have made the 
choice that the constraint should take precedence. 

We compared with the evolving constants approach and noted that the difference between 
that approach and the present one concerned the different ways in which equivalent classical 
expressions are turned into quantum operators. 

We showed that our class operators gave the correct and expected results when applied 
to standard non-relativistic quantum mechanics written in parametrized form. These results 
also showed how spacetime questions in non-relativistic quantum mechanics can be expressed 
in a fully spacetime form, involving an inner product defined on spacetime. 

We applied our formalism to some simple examples involving the free particle in one and 
two dimensions. These examples showed that the class operators could be easily calculated. 
Furthermore, there were some situations in which decoherence was possible for special states, 
without the need for an environment. However, the results were not always easy to interpret, 
and this is largely due to the fact that the classical limit of the quantum theory brings in 
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reflecting potentials. 

We briefly discussed the relativistic particle and considered the question of the probability 
of crossing a spacelike surface. We found that only states antisymmetric about the surface 
give decoherence and the crossing probability is then zero. This means that our formalism 
does not appear to reproduce earlier heuristic formulae for surface crossings. Furthermore, 
this result also emphasized that our approach is in fact a genuinely different proposal for 
the class operators compared to other approaches, and not a formalization of earlier more 
heuristic ideas. 

The formalism boiled down to particularly simple expressions for the case of systems of 
non-interacting harmonic oscillators and we exhibited some simple eigenstates which gave 
approximate decoherence and which had a clear semiclassical interpretation. 

A future publication will revisit the calculations of this paper but with the inclusion of 
an environment to produce decoherence js^j]. 
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